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the best method of simulta neously achieving simple structure and 
maximum similarity is the technique devised by Bloxora (1963). This 
technique simultaneously rotates different factor patterns to maximum 
similarity aod varimax simple structure and is applicable wher. the 
number of factors and variables are the same for each subpopulation. 
The technique is described and the computer algorithms are given. A 
numerical example compares the results of different factor analyses. 
If the limitations of the method are observed, Bloxon^s rotation 
technique has the potential for vide application in studies involving 
comparison of factor patterns. A small loss in simple structure 
through the use of Bloxom's method often pays off in factors which 
can be much more meaningfully compared. (CK) 
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Under some circumstances it is desirable to compare the 
factor patterns obtained from different factor analyses. When 
the number of factors and variables are the same for each sub- 
population, it is often possible to apply a rotation tectnlque 
developed by Blovora which simultaneously rotates different 
factor patxerns to maximum similarity and varlraax simple struc- 
ture. The technique is described, computer algorithms are out- 
lined, and a numerical example is given. 
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A review of a rotation to obtain maximum similarity and 
simple structure among factor patterns 

Carl Jensema 

University of Washington 
Introduction 

Factorial invariance is important in a variety of common research 
situations involving factor analysis. For example, an investigator 
may give the same questionnaire or test to several subpopulations and 
wish to know how stable the underlying factors arc across theso sub- 
populations, Alternatively, a researcher may give several similar 
questionnaires or tests to the same group and wish to determine vhether 
they all measure the same factors. 

In both of these two general situations or variations of them, 
the Investigator has a choice: He may rotate each of his factor pattern 
matrices to simple structure or he may rotate them to maximum similarity 
with each other. \ multitude of rotational techniques, both orthogonal 
and oblique, have been proposed for achieving simple structure transforms 
welch are unique for a given data set. Methods for achieving maximum 
similarity between factor patterns are somewhat less r.uraeious. They 
hk^ve developed along two general lines. The first is canonical correla- 
tions, which rotates correlational matrices to mutual maximum similarity. 
The second approach involves rotation to maximum similarity in relation 
to a target matrix. This later method is commonly known as Procrustes 



rotation. 



2 



As mentioned above, rotation techniques which are designed to 
achieve simple structure for a given data set tend to produce highly 
unique results. Comparison of different factor structure is usually 
difficult unless all factors are extremely similar. On the other hand, 
if rotation is ^.one to ach i eve maximum similarity between factor 
patterns, simple structure is almost always destroyed and meaningful 
interpretation is nearly impossible. Under certain conditions 
Procrustes rotation may produce results which are an exception to this 
general observations, but ever here, the target matrix is arbitrary and 
there is no way of ascertaining if it is the best possible one. 

It is highly desirable to have a rotational technique which will 
achieve both maximum simple structure and maximum similarity between 
different factor patterns. This has been recognized and several psycho- 
metricians have considered the problem. An article by Ahmavaraa (1954) 
provided some basic groundwork on the topic. Meredith (1964a) 
generalized Ahmavaraa' s work using Iawley's Selection theorem and 
proposed a rotational technique (Meredith, 1964b). Unfortunately, 
Meredith's solution has some rather severe limitations. It requires 
that the regression of the variates being factored on the variates de- 
fining the subpopulations be linear and homscedastic . The solution 
collapses if these conditions are not met. 

fcliff (1966) considered the problems of 1) rotating two solutions 
to maximum similarity and 2) rotating a solution orthogonally to a 
specified target matrix. He suggested a sequential procedure to obtain 
maximally similar simple structure solutions for two factor matrices. 

o 
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The problem with this treatment is that it may tend to be cumbersome 
and time consuming. 

To date, the best method of simultaneously achieving simple 
structure and maximum similarity is the technique of Bloxom (1968). 

He reversed Meredith 1 s approach and instead of assuming a constant 
factor pattern matrix while allowing the factor axore variance- 
covariance matrices to vary, he allowed the factor pattern matrices 
to vary across subpopulations while holding the factor .7" ere variance- 
co-variance matrix arbitrarily constant. This effectively circumverted 
the limitations in Meredith’s procedure while being more convenient 
and generally applicable than the treatment suggested by Cliff. 

However, Bloxom T s method appears to have escaped the general notice 
of psychologists. This paper will attempt to correct that situation 
by outlining the method and demonstrating its usefulness. 



where a is the factor loading matrix, H is the orthonormal transforma- 
tion matrix, and b is the rotated simple structure matrix. If there 
are ra subpopulation factor pattern matrices, each having the same number 
of variables and factors, the problem is to determine a transformation 
matrix, H, for each of the tn pattern matrices such that the resulting 
rotated matrices will be maximally similar and yet conform to the re- 
quirements of simple structure. 



General Method 



Rotation of factor patterns is usually described as: 



( 1 ) 
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The Varimax method of rotation, as described by Horst (1965), maximizes 



( 2 ) 
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where b is a matrix of the squared elements of the rotated loadings 

( l 

and \I - — J Is an indempotent matrix. Bloxom has generalized this 
criterion across ra loading matrices. His index of similarity between two 
factor loading matrices, ^b and ^b, is given by 
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as an expression of similarity between all m factor loading matrices and 
the extent to which the varimax criterion is maximized for each of the m 
loading matrices. By maximizing (4) with the constraint that 
(5) = jHjH 'a I, 



maximum similarity and maximum simple structure are achieved among the 
m rotated factor loading matrices. 
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Computat ion 



Initially, each of the k x k transformation matrices are set to 
identity matrices. The n x n inderapotent matrix, , is cal- 

culated by setting each diagonal element to - and each off 
diagonal element to — where n is the number of variables. After these 
initial * ^at ions are completed the iteration procedure is begun. Taking 
the unrotated loadings, ^a, of each of the ra matrices, the rotated 



matrices, ^b, are obtained from 

( 6 ) 



,0 - l* 1 H 



« I/ 2 ’ V / 2) l 



Squaring each element in each ^b matrix, we calculate 
ra m 

(7) * a ) )■ 

i-1 j-1 

where V is the inderapotent matrix, [l - . This is the criterion 

we seek to maximize. Iterations will be halted when \Jr stabilizes. 

Continuing the iterative procedure, equations (8) through (11) are 
looped through m times, once for each transformation matrix. Calculating 
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multiplication by the unrotated loading matrix, ^a, gives 

«) iV' • 

All roots and vectors of the minor product moment of this matrix are 



obtained, as indicated by 
( 10 ) 
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From this the new estimate of the transformation matrix is calculated as 



(11 



i H - U* *)(>' > cf v) • 



As previously indicated, equations (8) through (11) are carried out 
for each of m loading matrices. When this has been accomplished, 
equation (6) is returned to and a new iteration is begun. This process 
is continued until successive iterations indicate that has stabil- 
ized and final rotated matrices ^b are the best possible. 

Once the maximally similar varimax patterns have been found, it 
is usually desirable to have some measure of the degree of factorial 
similarity that has been achieved. The most common statistical index 
of similarity is the cosine between two factor vectors. The relation- 
ship between two factor pattern matrices, A and B, is given by 



(12) AT a B + E 

where T is a transformation matrix and E is a matrix of error terms. 

The T matrix may be viewed as an estimate of the cosines between factor 
pattern matrices A and B, and, as such, can be considered an index of 
factorial similarity, Schonemann (1966) has given a solution for T 
utilizing the constraints that 

(13) TT = I T » I 

and 

(14) tr (E / E) * min. 

Briefly, Schonemann’s T solution is calculated in the following manner: 

(15) S e A'B, 

(16) S'S « VD V t 

s 
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(17) S S = W D W ' 

s 

and 

(18) T = W V / . 

The W matrix of eigennectors in (17) must be reflected before (18) is 
computed. This may be done by obtaining the diagonal 

(19) = W'SV 

s 

1 X 

and examining D 2 for negative elements* If an element of D z is negative, 
s s 

reflect the corresponding column in W. When computed as described 
above, the T matrix has been found to be a useful index of factorial 
similarity for evaluating the output of Bloxora's rotation. 

Application 

Semantic differential data were used to demonstrate the usefulness 
of Bloxom's method in comparing the results of different factor analyses. 
Each of 176 subjects responded to the same 25 semantic differential 
bipolar word pairs for each of 3 concepts. These concepts were "Blind 
person", "Deaf person", and "Amputee". A principle components factor 
analysis was performed on each concept and three factors per concept 
were extracted. The factors for each concept were rotated with varimax 
and then with Bloxom's rotation. Table 1 gives the factor patterns 
obtained with each of these methods. Visual examination of the Bloxom 
rotation patterns readily demonstrates the close similarity between 
the three factors involved in each concept. This similarity is not 
nearly so apparent in the varimax factor patterns. 




8 



8 



Further evidence of similarity between factors is given by the 
cosines in Table 2. These provide a i^re precise measure of factorial 
similarity than simple visual examination. Cosines can take values 
ranging from +1.00 to -1.00. A cosine of +1.00 indicates perfect align- 
ment of factor vectors, .00 represents orthogonal factors, and -1.00 
means the poles of the factor vectors are reversed. From Table 2 it 
can be seen that the Bloxom method of rotation does an excellent Job 
of producing factor patterns which can be compared across concepts. 

As with any statistical technique, Bloxcm's rotation has limita- 
tions. Bloxom has pointed out that in some cases the method gives 
preferential treatment for facias having a large number of high load- 
ings. He suggested that this might be caused by a covariance "tug-of- 
war" between larger and smaller factors. If this happens, the factor 
patterns may tend to be variraax- covarimax for larger factors and 
varimax for smaller factors. Caution must be used in applying the 
method where factors within subpopulations are not of comparable 
magnitude . 

If the limitations of the method are observed, Bloxom 1 s rotation 
technique has the potential for wide application in studies involving 
comparison of factor patterns. A small loss in simple structure through 
the use of Bloxora*s method often pays off in factors which can be much 
more meaningfully compared. 
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Table 1 
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Factor Patterns Obtained with Two Rotation Methods 
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Cosines Between Factors 
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Program Setup Instructions 



A. Title Card — Any alphanumeriv message to be printed at beginning of output/ 

B. Problem Card 

Column Item 



1 - 7 
8-9 
10 - 11 



Punch "PROBLEM" 

Number of matrices (Max = 10) 

Number of variables in each Matrix 

(Max * 50) 



12 - 13 

14 - 16 
17 - 21 
22 - 23 



Number of factors in each Matrix 

(Max = 15) 

Maximum number of iterations to be allowed. 
Iteration tolerance (try 00100 or 00010) 
Number of input unit 



C. Input Format Card. Any F-type format 

D. Data — Matrices should be read in one at a time with factors arranged in 
columns and variables in rows 

This is a repeating program. Steps A through D can be repeated as often 
as desired. 



E. At the end of all problems put a card with "FINISH" punched in columns 1-6. 
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nno' 






i 



RUN ( S ) . 
SETCORE . 
LGO. 



-00 




2 



C 

3 



A 

20 

5 



PROGRAM MAXSlMt I NPUT .OUTPUT .PUNCH , T A^E 5= I NPUT » TAPE6- y UTPUT . TAP£7= 
♦PUNCH) 

DIMENSION A(10,!>0*15)»n(10»50,15) ,BV< 10 » 15 »50 > ,H (10» 15. lb) * 

*V( 50 .50) ,BB ( 50 .15 ) *AR( 15*15) »8AAE ( 30 * 1 5 ) .E < 15 1 , T I TLE ( 8 ) » FMT ( Q ) 
this program wap written by carl jensema in june.197o. 

IT TRANSFORMS FACTOR LOADING MATRICES OF DIFFERENT GROUPS TO 
MAXIMALLY SIMILAR orthogonal SIMPLE STRUCTURES. 
the method rs OUTLINED by bloxom in p.sychometrika 

(VOL. 33.NO. 2, JUNE* 196S). THE SIMPLE STRUCTURE CRITERION 
;s THE SAME AS FOR. VAR I MAX TRaNSFOKMAT I ONb. 

M=NUM8ER OF GROUPS (MAX=10( 



NV=N UMBER OF VARIABLES <MAX=50) 
NFsNUMBEK OF FACTOR* CMAX=15) 
A-ORK. I NAL FACTOR LOADING MATRICES 
H* TRANSFORM/ TION MATRICES 



8 = N E W FACTOR LOADING MATRICES 

READ IS • 1 ) TITLE. PROBLEM. M.NV.NF.MAX.TOL.NTtFMT 



FORM AT ( t A10/A7 *3 1 2 * I 3 » F5 • 5 ♦ 1 2/8A10 ) 

T F C TITLEIl ) .FO.l OH FINISH ) GO TO 10000 

IFIFMTIl I.EO.TOHFINISH ) GO TO 10000 
IF IPR03LEM.EQ.7HFI NISH ) GO TO 10000 
IF (PROBLEM. NE.7HPR0PLEM) GO TO 10 
WRI TE<6*2» TITLE.M.NV.NF.MAX.TOL.NT.FMT 
F0KMATIASH1 UNIVERSITY oF WASHINGTON BUREAU OF 
1//76H PROGRAM TU TRANSFORM FACTOR LOADINGS TU 
*MPLE STRUCTURES 2/1X.8A10 
2//19H NUMBER OF GROUPS = 17X.I3 
3//36H NUMBER OF VARIABLES IN EACH GROUP = 13 
4//36H NUMBER OF FACTORS IN EACH GROUP * 1 3 

5//36H MAXIMUM NUMBER OF ITERATIONS = 13 

6//22H ITERATION TOLERANCE = 1AX.F 9.5 
7//13H INPUT UNIT = ?HX.I3 
8//15H INPUT FORMAT = 10X.8A1G) 

READ IN ORIGINAL LOADINGS. 

DO 20 I — 1 . M 



TESTING 

MAXIMALLY SIMILAR SI 



WRITE (6. 3) I 

FORMAT ( 5 ( / J 1 28H ORIGINAL LOADINGS FOR GROUP 1 3 > 
DO 20 J= 1 >NV 



READ (NT.FMT ) ( A( I » J»KI .K=l .NF 1 
WRI T£l 6 * A ) U 1 ( A < ! . J > K ) » X= 1 . NF ) 

FORMAT ( //9H VAR I ABLE 1 3 » ( / 15F9.3)) 

CONTINUE 
WRI TEC6.5 I 

FORMAT ( 5 ( / ) * 32H COSINES BETWEEN INPUT MATRICES. ) 

CALL COS 1 NES ( A , BAAR • AB :E .TOL *rt , NV . NF » 

WRITE <6. A) 

FORMAT ( 5 ( / ) * 38H SUCCESSIVE ITERATION CRITERION VALUES > 
CALCULATE IDEMPOTENT MATRIX USED IN CRITERION FUNCTION. 

THE CRITERION IS THE SAME AS THE ONE USED IN VARIMAX TRANSFORMS. 
DO AO 1 = 1 .NV ' 

DO 30 J= 1 »NV 

Vf 1 .J) =0,0-1 ,0/NV 



14 



o 



:o 

40 

c 



• 50 
.60 



C 

c 

c 

70 



80 

c 

c 






100 

c 

7 

C 

*c 

r 



o 

ERLC 

iifflmaffaHaaaa 



CONTINUE 

V ( I , I ) = 1.0-1.0/NV 
CONTINUE 

SET TRANSFORMATION MATRIX H TO AN IDENTITY MATRIX, 

DO 60 1=1 »M 
DO 60 J=1 . NF 
DO 6 U K=1 >NF 
H ( I ,J,K)=0.0 
CONTINUE 
H( I » J » J ) = 1 • 0 
CONTINUE 
L00P=-1 
PSI=C.O 

BEGIN ITERATION LOOP. 

MULTIPLY ORIGINAL LOADINGS BY TRANSFORMATION MATRIX* 

THIS IS EOUATION 1 OF 8L0X0MS ARTICLE. 

L00P=L00P+1 
DO 80 1=1 *M 
DO 80 J= l »NV 
DO 80 K= 1 »NF 
B( I »J, 0=0.0 
DO 80 L= 1 »NF 

8( I » J,nJ+A( l* »L.iO 

CONTINUE 

BEGIN CALCULATION OF CRITERION FUNCTION PSI. 

THIS IS GIVEN IN EQUATION 5 OF 6L0X0MS ARTICLE. 

DO «0 I = 1 *M 
DO 90 J= i »NF 
DO 90 K=1»NV 
BV( I ,J.K) = C.O 
DO 90 L= 1 »NV 

6V<I»U.K)»8V< I .J.KI+I6I I .L ,U) **21 *V< L. K) 

CONTINUE 
OLDPSl =PS I 
PSI=C.O 
DO l.vl) I =1 *M 
DO 1.00 J = 1 » M 
DO 1 00 K=1 .NF 
DO 1 00 L= 1 >NV 

PSI =PS I +BV( I.K,L)*(B(J»L»K)*#2> 

CONTINUE 

PSI IS NOW THE CRITERION FUNCTION VALUE SPECIFIED BY EQUATION 5. 
WRI TEt 6 » 7 T PSI 
FORMAT! FI 5 .6 ) 

CHECK TO SEE IF ITERATION TOLERANCE HAS BEEN REACHED. 

IF( ABS(PSI-OLDPSIi .LE.TOLI GO TO 200 

CHECK TO SEE IF THE MAXIMUM NUMBER OK ITERATIONS HAS BEEN REACHED. 
IF( LOOP. GE. MAX) GO TO 200 

BEGIN INNER LOOP TO GET A NEW TRANSFORMATION MATRIX FOR EACH GROUP 

DO 190 11*1 »M 

DO 110 1 = 1 .NF 

DO i 10 J=l »NV 

8ft ( I » J ) =0 .0 

CONTINUE 

DO 1 20 1=1 »M 

DOUdLE=1.0 
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WR1 TE(6» 11 ) 

11 FORMAT { 5 ( / ) » 33H COSINES BETWEEN OUTPUT MATRICES. ) 

CALL COSlNES(BfBAAB>AB>E»TOL»M»NV.NF) 

GO TO 10 

10000 WRI TE ( 6 * 12) 

12 FORMAT ( 20 ( / ) , ] ?H END OF JOB.) 

STOP 

END 

SUBROUT 1 NECOSINES ( B »BAAB » AB »E » TOL »M , NV *NF ) 

C THIS SUBROUTINE CALCULATES ORTHOGONAL PROCRUSTES TRANSFORMATION 

C (COSINE) MATRICES. THE METHOD COMES FROM SCHONEMANN (PSYCHOMETRI KA » 

C MARCH, 1966). 

DIMENSION B 1 10, au* 15 ) » BAA8( 3^* IB) * ABBA! 30,15 ) * AB( 15* 15 ) »E( 1M 
DO 70 I 1 = 1 ,M 
DO 70 JJ= I I *M 

C OBTAIN MINOR PRODUCT OF TWO LOADING MATRICES. 

DO 1 0 1=1 ,NF 
DO 1 0 J= 1 »NF 
ABI I »J) =0.0 
DO TO K= 1 » NV 

AB ( I s J ) =AB ( I » J) +8 ( I I »K,I ) *R ( JJ * <» J ) 

10 CONTINUE 

c obtain minor and major products of ab. 

DO 20 I = 1 » NF 
DO 20 J= 1 » NF 
BAAd ( I ♦ J ) »0.f> 

A8BA ( I * J ) =0 .0 
DO 20 X* 1 *NF 

BAAdt I » J)=BAAR{ I » J ) +AB ( K » I )*ABl<* J> 

ABBA ( I * J ) = ABBA l I * J ) +AB ( I » K ) *AB( J » < ) 

20 CONTINUE 

C OBTAIN ROOTS AND VECTORS. 

CALL JACSIMIBAAB,E,TOL,NF) 

CALL JACSIM(ABBA.E«t0L,NF) 

C BEGIN REFLECTION OF ABBA VECTORS. 

C (SEE P • 8 » BOTTOM PARAGRAPH OF SCHONEMANN ARTICLE.) 

C MULTIPLY TO OBTAIN WSV. 

DO 30 1=1 »NF 
DO 30 J=1 *NF 
ABBA ( I »J >=0.0 
DO SO K=l»NF 
KNF=K+NF 

ABBA t I « J ) =ABB A ( I » J ) +ABBA. ( KNF * I )*AB(X*J) 

30 CONTINUE 

DO AC 1=1, NF 
DO AC J=1»NF 
AB ( I ,J)=0.0 
DO AO K=1»NF 
KNF=K+NF 

AB< I *J)=AB11 *J)+ABBA( I ,K)*8AAB(<NF» J) 

AO CONTINUE 

C ■ CHECX DIAGONAL UF WSV, REFLECT COLUMN IF DIAGONAL ELEMENT IS NEGATIVE, 
DO 50 I *1 . NF 

I F (ABC 1*1) . GE .0. ) GO TO 50 
DO 5 u J= 1 » NF 
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1F(!.EQ.II> DOUBLE “ 2.0 
DO T 20 J = 1 » N F 
DO 1 20 K=1 ,NV 
TEMP=E( 1 1 »K»J) 

DO i 20 L = 1 »NV 

B8U.J) = BB(<, J) + TEMP#VOv.L)*<B< I » L . J ) **2 I *OOUBLE 
120 CONTINUE 

.C WE NOW HAVE THE BETA MATRIX GIVEN IN EQUATION 18. 

C MULTIPLY THF BETA MATRIX BY THE ORIGINAL LOADINGS* 

c Then calculate the minor product moment of this new matrix, 

DO 1 '*0 1=1 .NF 
DO 1 40 J=1 »NF 
AB< I >J) = 0.0 

00 i AO Ks 1 , AiV 

AB( I »J)=AS< I »J)+A< II iK.I I*BBU»J) 

14D CONTINUE 

DO 1 50 1=1 »NF 
DO 1 50 J= 1 »NF 
BAAb ( I » J ) =0 .0 
DO 1 50 K=1 ,NF 

BAAiid *Jj = BAAR(I > J)+AB<<* I )*AB<K> J) 

150 CONTINUE 

C DO A COMPLETE FACTORIZATION OF MATRIX BAAB. 

CALL JACSIM(BAAB»E*TOL*NF) 

C EIGENVECTORS ARE COLUMNS OF BAAB. 

C EIGENVALUES ARE IN E. 

DO i 60 1=1 »NF 
E(I) =1 » O/SORT {£<!!) 

160 CONTINUE 

DO WO 1=1 .NF 

1 NF= I + NF 

DO i 70 J=1 .NF 
JNF = J=-NF 
BAABd . JI=0.0 
DO 1 70 K=1 »NF 

BAAol I » J ) =8AAR( I . J ) +BAAB 1 1 NF . K ) *E ( K ) *B AAb< JNF t<) 

17 0 CONTINUE 

DO 180 1=1 .NF 
DO 180 J=1 .NF 
HI II .1 . Jl =0.0 
DO 180 K=1»NF 

HI I I . I » J)=H( 1 1 » I . Jl+AB ( I «K)*UAAB< K. J ) 

180 CONTINUE _ 

C WE NOW HAVE NEW ESTIMATIONS OF THE TRANSFORMATION MATRICES. 

C END OF INNER LOOP 

19* CONTINUE 

' C END OF ITERATION LOOP 

GO TO 70 

-200 WRITEI6.8) LOOP 

8 FORMAT (// »}UH NUMBER OF ITERATIONS NEEDED = 15! 

. DO 210 1=1. M 

WRITEI6.9) I 

9 FORMAT (5 ( / ) »31H TRANSFORMED LOADINGS FOR GROUP 13) 

DO 210 J=1 .NV 

WRITEI6.4) J. (8(I»J«K)*KaltNF) 

210 CONTINUE 

O 
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JNF-J+NF 

A3 BA ( JNF > I ) r ABBA ( JNF » I > * ( -1 . ) 

CONTINUE 

END OF REFLECTION. 

obtain transformation matrix by multiplying vectors. 

DO 60 I =1 >NF 
INF= I+NF 
DO 60 J= 1 »NF 
JNF=J+NF 
ABI I »J)=0.0 
DO 60 K.= 1 sNF 

AQ ( I * J ) = AB ( I * J ) +ABBA ( I NF » K ) *BAAB ( JNF * K ) 

CONTINUE 

WR1 TE< 6* 1 ) II*JJ 

FORMAT! ///22H COSINES BETWEEN GR0UPI3*17H (ROWS) AND GROUPI3 
* 1 1H (COLUMNS). I 
DO 70 I = 1 »NF 

WRI TE(6*2 ) I » ( AB I i * J I » J=1 »NF ) 

FORMAT ( //7H FACTOR* 131 /15F9.3) ) 

CONTINUE 

RETURN 

END 

SUBROUTINE JACSIMt R»D»P*N) 

this subroutine calculates all roots and vectors of a matrix 

USING A JACOBI METHCN. 

DIMENSION RI30.15) »0( 15) 

Nl=N+) 

N]1=N-1 

N2=N*2 

DO 10 I=N1*N2 
DO 1 0 J=1 »N 
R ( I * J ) = 0 , 0 
CONTINUE 
DO 20 I = 1 » N 
NI-N+I 
R<NI»I)=J.O 

continue 

DO 00 L= 1 * 1 00 
DO 3 0 I = 1 ♦ N 
Dl 1 )*R< I »I ) 

CONTINUE 
DO 5 0 I = i *N1 1 
11= 1+1 

DO 50 J=I1»N 
DR=R 1 1 » I ) ~R I J * J ) 

AsSORT(DR**?+A.*R( I * J ) **? ) 

A=SQRT ( (A+DR)/(2.*A) ) 

B»S0RT( 1 ,-A«»? ) 

CsSIGNf 1. *R( 1 *J) ) 

DO A 0 K= 1 »N? 

U=RU.I ) *A*>C+R ( X * J I *8 

R(X.J)=-RIK>1 >*B*C+RI K* J)*A 

R(X»1)*U 

CONTINUE 

DO 50 K*1 *N 
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U=RU ,K)*A*C+R< U.K )*8 
R( J ,K)=-R< I *K ) *8*C+R ( J » K ) *A 
R ( I »K. ) =U 
50 CONTINUE 

DO 60 1 = 1. N 
DI I )=A8SID< I )-R< I. ID 
• 60 CONTINUE 

s=o < o 

00 70 1=1. N 

S=MAXl (S.D( I) ) 

70 CONTINUE 

DO 80 1=1. N 
D(II=R.'MI 
80 CONTINUE 

IFK-P) 100,100.90 
90 CONTINUE 

100 RETURN 

END 
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